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A simple formula for the infrared singular part of the integrand of one-loop QCD 
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We show that a well-known simple formula for the explicit infrared poles of one-loop QCD am- 
plitudes has a corresponding simple counterpart in unintegrated form. The unintegrated formula 
approximates the integrand of one-loop QCD amplitudes in all soft and collinear singular regions. 
It thus defines a local counter-term for the infrared singularities and can be used as an ingredient 
for the numerical calculation of one-loop amplitudes. 
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INTRODUCTION 



The calculation of ncxt-to-leading order (NLO) QCD 
corrections for multi-jet final states plays an important 
role for the experiments at the LHC. For any NLO cal- 
culation there are two parts to be calculated: the real 
and the virtual corrections. As far as the real corrections 
are concerned it has become a standard approach to first 
subtract out the infrared divergences |l|43| and to per- 
form the phase space integration of the subtracted real 
correction term numerically by Monte Carlo techniques. 
The subtraction terms can be integrated and added back. 
The integrated form has a remarkable simple structure 
and predicts the infrared poles of the renormalised one- 
loop amplitude. The one-loop amplitude contributes to 
the virtual part and has in general ultraviolet, soft and 
collinear divergences. It involves the integration over one 
unconstrained momentum (the loop momentum) in D di- 
mensions. It is tempting to follow a similar path also for 
the virtual part: First subtract out all singular terms 
of ultraviolet and infrared origin from the one-loop am- 
plitude and then perform the integration over the loop 
momentum together with the phase space integration in 
four dimensions by deforming the loop integration con- 
tour into the complex plane. The subtraction terms are 
then integrated analytically in D dimensions and added 
back. Such subtraction terms have been considered in 
[3| and issues related to the required contour deforma- 
tion have been studied in S]. The unintegrated in- 
frared subtraction terms in [J| are determined graph by 
graph. This limits the efficiency of the method. Inspired 
by recent work on the structure of infrared singularities 
of multi-loop amplitudes 043 show in this letter that 
the soft and collinear subtraction terms can be formu- 
lated at the level of amplitudes, without referring to in- 
dividual Feynman graphs. This is a significant simplifi- 
cation and opens the door to an efficient implementation 
based on recurrence relations [lo| . The required ultravi- 
olet subtraction terms have the form of propagator and 
vertex counter-terms and can be incorporated into the 
recurrence relations. A possible choice for the ultraviolet 



subtraction terms is given in ref. 

NOTATION 

Let us consider a one-loop QCD amplitude with n 
external partons. The full one-loop amplitude can be 
colour-decomposed into primitive amplitudes: 



(1) 



The colour structures are denoted by Cj while the prim- 
itive amplitudes are denoted by A^p ■ In the colour-fiow 
basis [11-1_3|] the colour structures are linear combina- 
tions of monomials in Kronecker (S^'s. Primitive am- 
plitudes are defined as a colour-stripped gauge-invariant 
set of Feynman diagrams with a fixed cyclic ordering of 
the external partons and a definite routing of the exter- 
nal fermions lines through the diagram 14. We have 
collected some information on primitive amplitudes in 
an appendix. It is convenient and sufficient to focus on 
primitive amplitudes. In the following we drop the sub- 
script j and focus on a single primitive one-loop ampli- 
tude which we denote by A^^' . Since the cyclic ordering 
of the external partons is fixed, there are only n different 
propagators occurring in the loop integral. We label the 
external momenta clockwise hy pi, p2, Pn and define 
<li ~ P1+P2 + ■■■+Pi, ki = k — qi. We can write the bare 
primitive one-loop amplitude in Feynman gauge as 
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bare 



d^'k (1) 
(27r)^^'"^'=' 

1 



kf — mf + iS 



(2) 



^iare ^^^'^ integrand of the bare one-loop amplitude. 
P{k) is a polynomial in the loop momentum k. The +i6- 
prescription instructs us to deform - if possible - the 
integration contour into the complex plane to avoid the 
poles at kf = rnf. If a deformation close to a pole is 
not possible, we say that the contour is pinched. In this 
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letter we restrict ourselves to non-exceptional external 
momenta. Then the divergences of the one-loop ampli- 
tude related to a pinched contour are either due to soft or 
collinear partons in the loop. These divergences are reg- 
ulated within dimensional regularisation by setting the 
number of space-time dimensions equal to D = A — 2e. 
A primitive amplitude which has soft or collinear diver- 
gences must have at least one loop propagator which cor- 
responds to a gluon. An amplitude which just consists of 
a closed fermion loop does not have any infrared diver- 
gences. We denote by Ig the set of indices i, for which 
the propagator i in the loop corresponds to a gluon. If 
we take the subset of diagrams which have the gluon 
loop propagator i and if we remove from each diagram of 
this subset the loop propagator i we obtain a set of tree 
diagrams. After removing multiple copies of identical di- 
agrams this set forms a Born partial amplitude which we 
denote by A -^^ . After integration, the soft and collinear 
poles of a primitive one-loop amplitude with massless 
partons are given by [H, [l5, 16 1 



^2pi ■ Pi 



(3) 



l,,2e 



We have multiplied the one-loop amplitude by S~ n 
where — (47r)^e~^^^ is the typical volume factor of 
dimensional regularisation, 7^ is Euler's constant and fi 
is the renormalisation scale. The constants Si are given 



by Sq 



Sn 



1 and Sg = 1/2, the index i of Si refers 



to the external particles. The generalisation of eq. ([3]) to 
massive partons is also known. 

Eq. ([3]) describes the infrared poles after integration 
over the loop momentum. In this letter we give a for- 
mula which approximates the integrand of eq. ([2]) in all 
infrared singular regions before the integration over the 
loop momentum. It thus defines a counter-term for the 
infrared singularities which is local in the loop momen- 
tum k and the external momenta pi. 



THE FORMULA 



part: 



(4) 



G 



(1) 

soft 



'ipi-Pi+l .(0) 

^1 I 



^ a 



^ WU. j ' ■ 

The Born partial amplitude J^i^ depends on the external 
momenta, but not on the loop momentum. The function 
guy ensures a regular behaviour of the collinear term in 
the ultraviolet region. A possible choice is 4] 



guv 



1 



2 V /^c kf /i^ 



(5) 



fic is an arbitrary scale. Integrating the soft and the 
collinear part we obtain 



d^k ^(1) _ as 



53 £.2 



{2Tr)D 47rr(l-£) 
2pi -Pi+i 



(6) 



S-^fi^ 



(27r)^ coi; 4^ p(j^ _ 



2 \ -e 



Equation (U) is the main result of this letter. This for- 
mula approximates the integrand of a primitive one-loop 
QCD amplitude in all soft and collinear limits. The ap- 
proximation is given by simple scalar two- and three- 
point functions, multiplied by a Born partial amplitude. 
One easily observes that the integrated form in eq. ([5]) 
agrees in the pole terms with eq. ([3]). 



PROOF OF THE FORMULA 

We first review briefly under which conditions infrared 
sin gula rities occur in an individual Feynman diagram [3, 
171 . Il8| : Soft singularities occur when a massless particle 
is exchanged between two on-shell particles. With the 
notation as in eq. ([2]) this corresponds to the case 



0, Pf 



2 2 

"^i-lJ Pi+1 



H+1- 



(7) 



In this section we present the formula for massless 
QCD. The extension to massive particles is discussed at 
the end of this letter. The infrared singular part can be 
written in unintegrated form as a soft and a collinear 



In that case the propagators (i — 1), i and (i -I- 1) are 
on-shell. The singularity comes from the integration re- 
gion k Qi. A collinear singularity occurs if a mass- 
less external on-shell particle is attached to two massless 



3 



propagators. This corresponds to 

= 0, rrii-i =0, rrii ^ 0. 



(8) 



In that case the propagators («— 1) and i are on-sheh. The 
singularity comes from the integration region k ~ qi—xpi, 
where x is a real variable between and 1. In order to 
proof eq. Q we consider now massless QCD amplitudes. 

The soft subtraction term is derived as follows: In the 
case where gluon i is soft, the corresponding propagator 
goes on-shell and we may replace in all Feynman dia- 
grams which have propagator i the metric tensor ij^i, of 
this propagator by a polarisation sum and gauge terms: 



2k^, ■ n 



(9) 



Here k\ denotes the on-shell limit of ki and d'^" denotes 
the sum over the physical polarisations: 

d^^{k,n) = Y,^l{k,n)e''_^{k,n) 

A 



2k-' 



(10) 



n'^ is a light-like reference vector. We note that self- 
energy diagrams are not singular in the soft limit, there- 
fore adding them to the loop diagrams will not change 
the soft limit. With the inclusion of the self-energy dia- 
grams and a corresponding replacement as in eq. (jH]) the 
contribution from the polarisation sum in eq. ([9]) makes 
up a tree-level partial amplitude, where two gluons with 
momenta k^ and —A;,- have been inserted between the ex- 
ternal legs i and i + 1. In the soft limit this tree- level 
partial amplitude is given by two eikonal factors times 
the tree-level partial amplitude without these two addi- 
tional gluons: 



P^+l ■ i-K) 



A 



(0) 



(11) 



In the soft limit we may replace 2pi-k\ by k'f_-^ and 2pi+i- 
{—k\) by kfj^^. Eq. PT|) then leads to the soft part of 
eq. (jl]). The terms with li^^ny and n^^k\'^ in eq. ^ vanish 
for the sum of all diagrams due to gauge invariance. 

The coUinear subtraction term is derived as follows: 
We have to consider diagrams, where two adjacent prop- 
agators in the loop go on-shell with a massless external 
leg in between. The cases where an external gluons splits 
into a ghost-antighost-pair or into a quark-antiquark-pair 
are in the collinear limit not singular enough to yield a di- 
vergence after integration. Therefore we are left with the 
case where an external quark splits into a quark-gluon 
pair and the case where an external gluon splits into two 
gluons. Let us first consider the q —> qg splitting. In 
Feynman gauge one can show that only the longitudi- 
nal polarisation of the gluon contributes to the collinear 



limit. The same holds true for the g gg splitting. In 
this case the collinear limit receives contributions when 
one of the two gluons in the loop carries a longitudinal 
polarisation (but not both). The external gluon has of 
course physical transverse polarisation. We can now use 
gauge-invariance to turn the sum of all diagrams with a 
collinear singularity in the propagators {i — 1) and i into 
the simple form of eq. (0]). The contraction of a longitu- 
dinal polarisation into a gauge-invariant set of diagrams 
yields zero. The set of collinear divergent diagrams forms 
an almost gauge-invariant set of diagrams. There is only 
one diagram missing, where the longitudinal polarised 
gluon couples directly to the other parton. This is a 
self-energy insertion on an external line, which by defi- 
nition is absent from the amputated one-loop amplitude. 
We can now turn the argument around and replace the 
sum of collinear singular diagrams by the negative of the 
self-energy insertion on the external line. The self-energy 
insertions on the external lines introduce a spurious 1 /pf- 
singularity. In order to calculate the singular part of the 
self-energies we regulate this spurious singularity by al- 
lowing p^ slightly off-shell, but keeping fci_i and ki on- 
shell and imposing momentum conservation. We can use 
the same parametrisation as in the real emission case: 



xp + fcj_ — 



-fci = (1 - x)p - fc^ - 



X (2p ■ n) ' 
kl 



(1 - x) {2p ■ n) 



(12) 



with p^ 



and 2p- k± = 2n- k± =0. The singular 



part of the self-energies with one longitudinal polarised 
gluon is proportional to 
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■plong 

9^99 ~ 




2 
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2ki, 


-1 


ki 


\ X 1 



2^i>, (13) 
+ 2 ) d^''(p, n). 



The terms with 21 x and 2/(1 — x) correspond to soft 
singularities and have already been subtracted out with 
the soft subtraction term G^^^j. In the collinear limit we 
therefore just have to subtract out the terms, which are 
non-singular in the soft limit. These terms are indepen- 
dent of X and lead to the collinear part of eq. (|4]). 

A few remarks are in order: Eq. ([3]) is usually stated 
for the renormalised one-loop amplitude. 



^(1) = - 



47r r(l - e) 



E 

pairs (jj) 



2T,T, 



-2k • p, 



I -1-1 



E 



(14) 



The colour charge operators are given by 

^j&ac -fpj. final state particles. 



T,- = -T- and Tg 
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The constants 7^ are given by = = 3Cf/2 and 
7g = /3o/2 with /3o = 11Ca/3 - 4TijiV//3. The renor- 
mahsed one-loop amplitude with Ug gluons, Ug quarks 
and rig antiquarks is related to the bare amplitude by 



(15) 



xAbare {pi, . ■ . , Pn, Zg S ^ V^^"s) • 



Zg is coupling renormalisation constant 



(16) 



Z2 is the quark field renormalisation constant and Z3 is 
the gluon field renormalisation constant. The LSZ re- 
duction formula instructs us to take for the field renor- 
malisation constants the residue of the propagators at 
the pole. In dimensional regularisation this residue is 1 
for massless particles and therefore the field renormalisa- 
tion constants are often omitted from eq. (jl5p . However 
■2^2 = ^3 — 1 is due to a cancellation between ultraviolet 
and infrared divergences. In Feynman gauge we have 



Z2 
Z3 



= 1 



— 0_F 

as 



— - — ]+0{al), (17) 



1 



1 + ^(2C, 



£UV 



1 



1 



Here we indicated explicitly the origin of the 1/e-poles. 
Using Z2 = = 1 mixes therefore poles of ultravi- 
olet and infrared origin. There are different variants 
of dimensional regularisation: conventional dimensional 
regularisation, the 't Hooft-Veltman scheme and four- 
dimensional schemes. The soft and collinear subtrac- 
tion terms in eq. (HI are independent of the variant of 
dimensional regularisation. The infrared poles obtained 
from the phase space integration of the real emission am- 
plitude depend on the other hand on the scheme of di- 
mensional regularisation. Again the solution is due to 
the field renormalisation constants which mix scheme- 
dependent terms of ultraviolet and infrared origin [l9| . 
The scheme-dependence of the bare one-loop amplitude 
is entirely of ultraviolet origin. 



GENERALISATION TO MASSIVE PARTONS 

In this section we present the generalisation to mas- 
sive QCD. This is in particular relevant to top quark 
physics. There are only a few modifications necessary 
with respect to the massless case. The modification for 
the unintegrated soft subtraction term is straightforward: 



'-'soft - 

—AnaA 



4pi -Pi+i 



-A 



(18) 

(0) 



As before, the sum runs over all particles in the loop 
which are gluons. Integrating the soft subtraction term 
we have to distinguish whether the masses m^-i and 
TOj+i are zero or not. The result can be written as 



1 ..2£ 



E 



{2tt)d 47rr(l-e)^ 
X C {{p, + p,+if,mj_i,mj+i, f^') Af^ + 0{e) 



(19) 



where the function C (s, rri^ , , n^) is given for the four 
different cases by 18|, |20| 



,2 ,,2\ _ /-</'„ ™2 n 



(20) 



C(s,m^O,M') = ( -J 



1 2, 
^ + - m . 

e e \ m'^ 



m 



-2Uo 



-21n(l -x^)+ln 



2x{s — ra\ — m^) 
77117712(1 — x^) 

77lim2 



ln(a;) - — +U2{x^) 
b 



ilnMl^ 

2 \m2 



L12 1 - X — + L12 1 - X — 

7712 / V "^1 



with 



1-X 



X 



1 



4mi77l2 



s — (7711 — 7712) 



(21) 



The modification for the collinear subtraction term is 
even simpler: There is no collinear singularity if an ex- 
ternal quark or antiquark is massive. It suffices therefore 
to define Sq — Sq — for a massive quark or antiquark. 

CONCLUSIONS 

In this letter we have shown that the infrared singu- 
lar part of the integrand of a primitive one-loop QCD 
amplitude is given by simple scalar two- and three-point 
functions, multiplied by a Born partial amplitude. The 
immediate application is to use this form as a subtraction 
term for the numerical integration of a one-loop ampli- 
tude. In view of refs. we do not exclude the possi- 
bility that this method generalises to higher loops. 

Primitive amplitudes 

In this appendix we include a brief summary on prim- 
itive one-loop amplitudes. In order to construct a prim- 
itive one-loop amplitude one starts to draw all possible 
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FIG. 1: Examples of diagrams: Diagrams (a) and (b) con- 
tribute to the left-moving primitive amplitude, while diagram 
(c) contributes to the right-moving amplitude. 



A diagram is translated to a formula with the help 
of the following Feynman rules: The propagators for a 
quark, gluon and ghost particle are 



fc2 



(22) 



respectively. The Feynman rules for the vertices are listed 
in fig. m 

In the main part of the paper we made use of the Ward 
identity for cyclic ordered amplitudes. The Ward iden- 
tity states if we replace in the amplitude the polarisation 
vector of one external gluon by its momentum then we 
obtain zero. A proof can be found in the textbook by 
Peskin and Schroeder 21 1. 




= [g"" {k^ - kf) + {kl - fc^ ) -h g'" (tr - t;)] 



= i[2g''V-g'"'g"'-g'"-g-'] 



FIG. 2: Colour-ordered Feynman rules. 



planar one-loop diagrams with a fixed cyclic order of the 
external legs, subject to the constraint that each fermion 
line is either only left-moving or only right-moving. We 
call a fermion line "left-moving" , if following the arrow of 
the fermion line the loop is to the right. In an amplitude 
with r external quark-antiquark pairs p of these pairs can 
be left-moving The remaining (r—p) pairs are then right- 
moving. As an example we consider primitive qggq- 
aniplitudes with the cyclic order (Ig, 2(7, 3g, 4,). Fig. [T] 
shows some diagrams as examples. Diagrams (a) and (b) 
contribute to the primitive amplitude, where the quark 
line is left-moving, diagram (c) however contributed to 
the primitive amplitude, where the quark line is right- 
moving. We would like to point out that diagram (b) has 
to be included. In an analytic calculation this diagram 
is often discarded, because it yields zero within dimen- 
sional regularisation. However this zero is obtained from 
a cancellation between an ultraviolet divergence and a 
collinear divergence. In a numerical calculation we have 
to keep this diagram in order not to spoil the local struc- 
ture. 
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